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A NJL-based study of the QCD critical line 
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We employ a 3 flavor NJL model to stress some general remarks about the QCD 
critical line. The dependence of the critical curve on fj, q = (fi u + Hd)j1 and /// = 

O ; 

(fj, u — Hd)/2 is discussed. The quark masses are varied to confirm that, in agreement 

(N . 

(3JT)' with universality arguments, the order of transition depends on the number of active 

flavors Nj. The slope of the critical curve vs. chemical potential is studied as 
a function of Nj. We compare our results with those recently obtained in lattice 
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simulations to establish a comparison among different models. 
PACS numbers: ll.10.Wx, 12.38.-t, 25.75.Nq 



INTRODUCTION 



Ph. In recent years, the study of the QCD phase diagram by means of numerical lattice 



simulations has improved considerably. In particular, thepresence of a critical ending point, 
first discovered within microscopical effective models DDD , appears to be a solid feature 
, althought its exact location along the critical curve is still controversial. Moreover, 
the incoming realization of LHC enhances the interest of the scientific community in this 
kind of problem. 

The main problem regarding the lattice study of QCD phase diagram is related to the 
so called sign problem; the fermionic determinant is not positive definite at finite baryon 
chemical potential, and therefore, to avoid this unwelcome feature, some suitable tricks 
are needed. The most commonly used between them are the study of QCD at imaginary 
chemical potential 0, Q, the reweighting procedure jj] and Taylor expansion in fi/T 
Q|. A possible way of avoiding the sign problem is to consider QCD at finite isospin 
chemical potential [ir. in this case, the fermion determinant is real and positive definite, 
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121 ]. Futhermore, the regime of finite 



and standard Monte Carlo simulations are allowed 
/i/ can be also studied within a class of effective models, and the results can be compared 
with those obtained on the lattice in order to check the consistency of different approaches. 
Two properties are remarkable and are of interest at /// 7^ 0: pion condensation and the 
splitting of critical curves related with light flavors. The former has been the subject of 
several studies within different models, starting from effective Lagrangians ^| , random 



matrices 



NJL 
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rj and ladder-QCD [lj and so apparently a model-independent 
phenomenon. However, pion (and kaon condensation would not be accessible from 
Heavy Ion experiments, but could regard the physics of compact stars. The latter question 
is slightly more controversial: the role played by instantons seems to be crucial to determine 
whether values of [ii obtainable in experiments can produce the separation of critical lines 

Q- 

The aim of this paper is to offer an overview of results concerning the behaviour of the critical 
line, obtained in the NJL model and directly comparable with recent lattice analyses. In its 
strong simplicity, the NJL model recovers the basic structure of non perturbative dynamics 
ruling the problem. Therefore, it can be trusted as a good toy model for the study of QCD 
phase diagram. 



II. THE MODEL 

Let us now consider the Lagrangian of the NJL model with three flavors u, d, s, with 
current masses m u = m d = m and m s and chemical potentials fi U i^di^s respectively 



where 



C — Cq + C m + + £4 + Cq 



a=0 

+ K [detail + 75)^ + det^(l - 75)*] 
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ftA is the current quark mass matrix which is taken diagonal and A, is the matrix of the 

[2 

quark chemical potentials. As usual Ao = y - I and A a , a = 1,...,8 are the Gell-Mann 
matrices. 

The 't Hooft determinant term, for the three flavors case, corresponds to a six fermion 
interaction. By working at the mean field level, the six fermion term can be recast into 
an effective four-fermion one. In such a way the Lagrangian (J2J) reduces to the usual NJL 
Lagrangian, apart from a redefinition of the four-fermion coupling constant G into a new 
set of effective ones, taking into account the flavor mixing arising from the 't Hooft term 
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22] 



Therefore, because we are dealing with four fermion interactions only, we can calculate 
the effective potential by using the standard technique to introduce bosonic (collective) vari- 
ables through the Hubbard- Stratonovich transformation and by integrating out the fermion 
fields in the generating functional. 

If we limit ourselves to consider the three scalar condensates, and the pseudoscalar conden- 
sate in the light quark sector only, the one-loop effective potential we get is: 

V = ^q(xI + X 2 d + xl) - A3 J^ Xu Xd Xs + ^(1 - ^^r)(pld) + ^log ( 3 ) 



where 



n=+oo 



log= p J (2vr) 



( h u -7 5 A Pll(i (l-^) 0^ 

75 A p ud (1 - ^) h d 



V 



8G 2 

h 



s 



(4) 



h f = (iuj n + ///)7q - p-7 - Mf 



In eq. (@J) Tr means trace over Dirac, flavor and color indices and u n = (2n + 1)tt//3 are the 
Matsubara frequencies. The dimensionless fields Xf an d Pud are connected to the scalar and 
pseudoscalar condensates respectively by the following relations 
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x , = - 4G <*'*'> 



Pud — — 2G 



A 

(5) 

(uj 5 d - dj 5 u) 
A 



and are variationally determined at the absolute minimum of the effective potential. The 
constituent quark masses are 

M t = m l + A Xl -A 2 ^^ (i^j^k) (6) 

Since the model is non renormalizable, we have to introduce the hard cut-off A on the 
three-momentum. 

Through this letter, when we consider the physical case (with realistic values of meson 
masses and decay consant), we assume for the parameters the same values as in ref. (3] 

A = 631.4 MeV; G A 2 = 3.67; K A 5 = -9.29; (7) 



rh = m » + md = 5.5 MeV; m s = 135.7 MeV (8) 

To investigate regimes different than the real one, with a varying number of Nf massless 
flavors, we will treat quark masses as free parameters, by keeping coupling and cut-off scale 
fixed as in eq. (J7J). 

In the following, it will turn out to be more convenient to introduce the following linear 
combinations of chemical potentials: 

fi q = (pu + Pd)/2] fa = (fa - JJ>d) (9) 



The quark chemical potential n q is just one third of the baryon chemical potential ii q = 

yses previously performed within the chiral Lagrangian approach 

studies, we expect a superfluid phase with 
condensed pions when the isospin chemical potential Hi exceeds a critical value (pf = 
mn/2 at /i g = T = 0). 



a 3. Following the analyses previously perfor 
Q, and NJL 17, 23 1 and ladder-QCD Q 
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III. BEHAVIOUR OF CRITICAL LINES 

A. Critical temperature dependence on baryon/isospin chemical potentials 

The aim of this section is to shed some light on the physics of QCD at finite baryon 
chemical potential (fi q = /ig/3), by comparing the physics at fi q ^ and \ii = with that 
at \i q = and fij ^ 0. In fact, the latter case can be studied on the lattice by means of 
standard importance sampling techniques. The connection of these two regimes could give a 
deeper understanding of the sign problem in the fermion determinant, and provide us with 
some procedure to check present simulations and possibly improve numerical algorithms. 
Rigorous QCD inequalities at non zero chemical potential have been proposed to try to 



resolve this question in |24| : the enigma why, at T = 0, there exists a critical value for 
chemical potentials below which the system lies in its ground state has been called "silver 
blaze problem". A possible solution of this problem, by using 1/N C expansion, has been 
proposed in j3, Q • A recent study of phase quenched QCD (a theory where the absolute 
value of the fermion determinant is taken) has been performed in j^ . 
We will consider here the dependence of the critical temperatures T c {ji q ) = T c (fi q ; fii = 0) 
and T c (/x/) = T c (/xj; fi g = 0), for low chemical potentials, obtained by a mean-field analysis 
of the NJL model. Obviously, mesons and baryons (and di-quarks) carry different spin and 
charges, and their properties depend differently on fi q and ///•; for these reasons, one could 
expect the two curves T c (/i q ) and T c (/ii), when starting by the same value at \i q = fii = 0, to 
be different, at least in the regime where bound states heavily influence the thermodynamics 
of the system. On the other hand, when the free-energy is mainly ruled by the constituent 
quarks, there is no reason to expect a dependence of the critical curve on the signs of chem- 
ical potentials (which will fix the sign of the total charges associated with the system). 
Here at the mean-field level, the effect of bound states is considered when we admit the 
formation of a pion condensate. Actually, in agreement with chiral models analyses, in the 
NJL mode 
computed 

of the isospin chemical potential. For this reason, as /i/ is higher than some critical value 
(m T /2 at T = 0), a pion condensate starts to form; a similar effect happens when \i q is 



the pion effective mass dependence on chemical potentials can be analitically 
2^|; the result is that the charged pions chemical potential is exactly the double 



higher than the critical value for di-quark condensation, which is expected to occur at val- 
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ues fi q > 300 -7- 400 MeV (of course, before di-quark condensation, for ji q ~ m N /3 and low 
temperatures there should be the liquid-gas transition for the nucleons) . Since in this paper 
we are interested mainly in the regime of relatively small chemical potentials (lower than 
~ 200 MeV) we will neglect the latter possibility. 

In fact, when the pion condensate is zero, the mean field effective potential is symmetric 
under [i u — > —fi u , fid —> ~^d] this implies that ji q <-> ±/i/ is a symmetry of the problem. 
Therefore, for zero p, the two curves T c (/i q ) and T c (/ij) have the same analytical dependence. 
In Fig. [H |21 the phase diagrams in (p q ,T), (/i/,T) spaces are shown; starting from the com- 
mon value T = 201 MeV, corresponding to p q = /if = 0, the cross-over curves coincide up to 
the value of about 150 MeV for both chemical potentials. For higher [ij (and temperatures 
lower than ~ 200 MeV) we are in the condensed pions phase; in agreement with 23] this 
regime will persist until fii ~ 860 MeV, before the saturation regime takes place. On the 
other hand, if we follow the cross-over line T c (/i 9 ), we will find, as expected, a critical ending 
point for fi q = 330 MeV, T = 42 MeV, and a line of first order transitions for higher \i q . 
Summarizing, at this level of calculation (mean field), until isospin chemical potential is 
lower than the critical value for pion condensation the curves T c (n q ) and T c (fij) have the 
same analytical expression T c (n q ) = T c (fij). This could be interesting in the attempt to 
extend lattice results from /i/ ^ and \x q = to p q ^ and /i/ = (at least in the region 
of low chemical potentials). Our conclusions appear to be in agreement with the authors of 
ref. |3] (both from the lattice and from a hadron renonances gas model) . 



It is generally assumed from universality arguments 29] and lattice analyses, [3flj . that 
the order of the phase transition by increasing temperature at zero density can change as the 
quark mass values are varied. If in the realistic case (with physical quark masses) the zero 
density transition is expected to be a cross-over, lattice analyses seem to show a first order 
transition when the three light quark masses are small enough. In particular, by taking 
m u = rrid = 0, there should be a critical value for m s below which the transition turns into 




B. Order of the transition by varying m s 




7 



250 
T (MeV) 

200 




100 150 200 250 300 350 

IH (MeV) 



FIG. 1: Phase diagram in the plane (fx q ,T), for the physical 2+1 case; fij and fj, s are set to zero. 
Dashed/solid lines indicate cross- over/ first order transitions; consequently, the dot in the picture 
labels the critical ending point. 



a discontinuous one: different lattice approaches find wR to 



)e half of the physical value of 



3 



22j . and we take 



the strange quark mass 

31] 

or rrig ~ 5 -T- 10 m u ^ (physical) 
To study this aspect in the NJL model, we start from the parameters fit of 
the quark masses as free parameters; namely, we take the four and six fermion couplings 
fixed so as to reproduce the phenomenology of the realistical physical situation. This way of 
proceeding could seem as rather arbitrary, but for any value of the masses considered here, 
we have verified that the output parameters have reasonable values (critical temperatures 
130 -i- 200 MeV, light quark scalar condensates (— 250 -: — 240 MeV) 3 , constituent light quark 
masses 250 -j- 350 MeV). Of course our results will be strongly model dependent, both for 
the choice of a specific set of parameters and of a particular model in itself. For instance, the 
NJL model provides an estimate of the position of the critical point at lower temperatures 
and higher chemical potentials with respect to those obtained in ladder-QCD and the ones 
from recent lattice simulation j^. 

In Fig. 01 we show the behaviour of the light quark scalar condensates vs. temperature, in 
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FIG. 2: Phase diagram in the plane (pi, T), for the physical 2+1 case; ji q and fj, s are set to zero. The 
dashed line indicates a cross-over transition for the scalar condensates, whereas dotted line stands 
for genuine second order transition for the pion condensate. In the regions of a non-vanishing 
pion condensate, the discontinuous behaviour of the scalar condensate turns into a continuous one, 
therefore the critical ending point is not present in this case. 

the m u = rrid = limit, and for vanishing chemical potentials too. In the upper picture, 
m s is taken to be zero, and for a temperature of about 130 MeV there is a sharp first 
order transition. As we increase m s the discontinuity of the scalar condensates reduces (in 
the lower picture the case m s = 8 MeV is shown), and when m s exceeds the critical value 
nig = 10 MeV, the zero temperature transition turns into a genuine second order transition. 
The value we get for is in any case smaller than the one from lattice predictions. 
In Fig. H we plot the phase diagram in the (fi q ,m s ) space, by taking m u = rrid fixed to 
zero, and fii = fi s = 0. The label //// indicates, for every couple (fi q ,m s ), whether, by 
increasing the temperature starting from zero, the transition is a first or a second order one. 
Actually, up to m s < = 10 MeV, we have first order transition for every value of \x q , 
and consequently there is no critical point; for m s slight above the critical value, the critical 
point locates at a value of fi q ~ 200 MeV. The critical value for fi q grows together with 
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m s until the strange quark decouples from the two light quarks and the critical point ii q 
coordinate is independent on m s , and lies at \i q ~ 300MeV. 

Finally we have studied whether, in agreement with lattice analyses, a first order transition 
persists when we consider a non zero but small m u = m^. this does not happen in the NJL 
model with our choice of parameters, for any value of m s . A recent work based on the linear 
sigma model had found the critical value for m u = = m s to be m crit = 40 ± 20 MeV • 
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FIG. 3: Behaviour of the light quark scalar condensates as a function of temperature at zero 
chemical potentials and m u = rrid = ; for m s = (upper picture) and m s = 8 MeV (lower 
picture). In the upper picture, X = Xu = Xd = Xs, in the lower X = Xu = Xd- When the strange 
quark mass exceeds the critical value = 10 MeV, the discontinuous behaviour turns into a 
genuine phase transition. 
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FIG. 4: Phase diagram in the (fj, q ,m s ) space, at m u = m<f = and fj,i = fj, s = 0. For every 
region in the diagram, the label I /II means whether, by increasing temperature and starting from 
T = 0, the transition is of first or second order. Actually, the line in the diagram separating the 
two different regions follows the critical point by varying m s . 



Even though lattice analyses at finite density still present ambiguities in their different 
approaches, some general features about the critical line appear to be rather solid. In 
particular, if To is the critical temperature for zero chemical potentials, the dependence of 
(T/T ) as a function of (fJ>/T ) should be parabolic, (at least in the regime (/x/T ) < 1), i.e. 
of the form (T/T ) = 1 — a(fi/T ) 2 . Secondly, the a coefficient should depend on the number 
of flavors Nf, increasing with Nf] in fact, the curves relative to Nf = 2, 2 + 1, 3 should be 
very close to each other and the one relative to Nf = 4 should be steeper. 
It is clear that a dependence of the a coefficient on Nf must be related in any case with 



C. 



Critical lines as a function of N 



f 
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a coupling between the flavors; otherwise, the effective potential would become a sum of 
single flavors contributions, and the critical temperature would not depend on Nf. This 
fact can give us an idea about the strength of the coupling between the flavors at the phase 
transition, and of its possible reduction with temperature. We will check the issue relative 
to the cases 2, 2 + 1, 3 in the framework of the NJL model. 

In this context, the aforementioned situations labeled by 2,3 are those with 2,3 massless 
flavors (in the case 2, m s is set to 5 GeV to decouple the strange quark); 2+1 is the physical 
case with realistic values of quark masses. In the following, // will indicate a common value 
for the chemical potential equal for all the active flavors. 



T/T 




FIG. 5: Plot of the three cases 2,2 + 1,3 in dimensionless units (ii/Tq,T/Tq), for h/Tq < 1. 
The curves relative to 2 and 2 + 1 (the dashed one) almost overlap. In disagreement with lattice 
simulations, the curve relative to 3 stands slight above the others. 

In Fig. we show the phase diagrams relative to the cases 2 + 1,2, 3. Obviously, in 
case 2, we are in the situation where m s < and we have only first order transitions. 
As we are varying the quark masses to consider different situations, Tq has a large range of 
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variation (from 130 to 200 MeV); it is impressive that, as we plot the phase diagrams in 
dimensionless units, those large differences cancel out almost completely. This is the most 
striking evidence that our approach, in its simplicity, has some validity. 

When showing the three cases on the same diagram, we can clearly observe the surprising 
overlap of the results. The agreement between 2 and 2 + 1 cases (apart from numerical 
instabilities related with 2 + 1 case) is remarkable, and good between 2 and 3 cases; the 
slope of case 3, with respect to the lattice previsions, is smaller than that of 2, but the 
difference is slight. 

The results we get for a are the following: 

2 a = 0.1995 + 0.025 

2 + 1 a = 0.2496 ± 0.0623 (10) 

3 a = 0.1614 + 0.011 

The worse precision we obtain for the 2 + 1 case depends on the larger error in the 
determination of the cross-over curve. 

For completeness, we quote the results for a from lattice analyses: 



2 


a 


= 0.0507 ± 0.0034; //„ = p: d = // [6] 




2 


a 


= 0.0504 ± 0.0036; fi u = fi d = fi [8] 




2 


a 


= 0.07 + 0.03; l i u = l i d = l i [10] 




2 + 1 


a 


= 0.0288 ± 0.0009; fi u = fi d = fi s = fj, 


= /W3 [4] 


3 


a 


= 0.0610 ± 0.0009; /x u = fi d = fi s = \i 


= /W3 [8] 


3 


a 


= 0.114 + 0.046; fx s = [10] 




4 


a 


= 0.099; fi f = fi = ^b/3 [6] 





(11) 
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We find that our predictions for a are bigger than those obtained by lattice approaches, 
apart from for Nf = 3; in that case, the results are comparable. A recent study based 
on a hadron resonance gas model |3| give the result a = 0.17 ± 0.01 and this value is in a 
good agreement with our results. A study within a chiral quark model give for a a value of 
about 0.1, extracted from Fig. 2 of 1 3 311 . 

n 

We have also studied the behaviour of the critical curve in ladder-QCD |3J (in its version 
]); in this model there is no coupling between flavors, therefore it is independent on 
Nf. We have found that the critical curve is flatter than that of NJL model, namely the 
coefficient a is much smaller, and hence closer to lattice predictions: a = 0.0797 ± 0.0056. 
We have attributed this feature to the lack of coupling between flavors in the model. On the 
other hand, since the value for the a coefficient we find in the NJL model is slightly higher 
than the value found in ref. (28] and sensibly higher than the values obtained from other 
lattice analyses, we can argue that introducing an effective reduction of the K coupling with 
temperature, as considered in ref. j3, S]> a would decrease at the meantime. Therefore, 
by considering the following temperature dependence of the 't Hooft term 

K{T) = K exp(-T/T 1 ) 2 (12) 

we have verified that reducing T\ the critical temperature at [A — is also reduced; in 
this way, the curve gets flatter. By taking two different values for T\, and considering for 
simplicity the case 2, we find the value a = 0.1186 ± 0.0061 in the case T\ = 160 MeV, and 
a = 0.1084 ± 0.0035 in the case T\ = 100 MeV. In this way, the agreement with lattice is 
considerably improved; this can be considered an indirect proof of U (1)^ effective restoration 
with temperature. 

For the sake of completeness, we have also studied the two flavor NJL model with the 
't Hooft determinant. The model is not completely equivalent to the three flavor model in 
the limit of infinite m s , since in the latter case a strange quark loop gives a contribution, 
proportional to the strange quark condensate, to the light quark constituent mass. In any 
case we do not expect a dramatic change in our results, with respect to previous expectations, 
for the behaviour of the critical curve. 

We consider in this case the following expression for the interaction part of the Lagrangian 
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(the U(1)a breaking determinant term can be rewritten as a four fermion interaction): 



dm = G 1 [{qq) 2 + {qTq) 2 + {qt l5 q) 2 + {qt l5 Tq) 2 }+G 2 [{qqf (13) 
with 

Gi = (l-/3)G ; G 2 = (3G (14) 

The j3 coefficient tells us how hard the flavor mixing is; it is maximal for G\ — 0, namely 
for (3 = 1. 

For the choice of the parameters G\ and G 2 we follow the approach proposed in ref. 
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2J. In ref. 



22j the authors study the original two flavor Lagrangian, proposed by 
Nambu and Jona-Lasinio, with G\ = G 2 and therefore (3 = 0.5. The value we find in this 
case for a is very similar with the result we obtained in the SU (3) case in the limit m s — » oo: 

a = 0.2107 ±0.0214 (15) 



The authors of ref. [2Cj take (3 as a free parameter instead. Here we furthermore 
consider the possible dependence of G 2 coefficient on the temperature, G 2 = G 2 (T = 
0) exp(-(T/T ] 



\2\ 



If we take G 2 independent on the temperature (namely with T\ = oo), the value we find 
for the a coefficient does not change by varying (3: 

a = 0.2142 ±0.0259 (16) 

again in agreement with previous analyses. On the other hand, if we admit a restoration of 
U(1)a symmetry with temperature, we find a dependence on the (3 coefficient. 
For (3 = 0.2 we have 

a = 0.1484 ± 0.009 for T x = 160 MeV (17) 

and 

a = 0.1196 ± 0.0161 for T x = 100 MeV. (18) 

For (3 = 0.3 we have 

a = 0.1024 ± 0.0042 for T x = 160 MeV (19) 
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and 

a = 0.08101 ± 0.007 for T\ = 100 MeV. (20) 

However, according to Shuryak 36J it is very unlikely that restoration of U{1)a can occur 
before chiral symmetry restoration; therefore, the value T\ = 100 MeV should not be taken 
too seriously. In any case, it appears clear that restoration of U(1)a symmetry can strongly 
influence the behaviour of the critical curve. 



IV. CONCLUSIONS 



In this paper we have studied some general features of the QCD critical line in the 
framework of a NJL model. In section one, we have compared the physics at fi q ^ with 
that at hi ^ 0. In section two, we have varied quark masses to show that the order of finite 
temperature transition changes if we consider small enough masses. In section three, we 
have studied the dependence of the slope of the critical curve T c vs. fi on Nf. In recent 
times these questions have received much attention from the lattice community, due to the 
great improvement of finite chemical potential algorithms of simulation. 
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